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Through Quantum Monte Carlo simulation, we study the biquadratic-interaction model with the 
SU(2) symmetry in two and three dimensions. The zero-temperature phase diagrams for the two 
cases are identical and exhibit an intermediate phase characterized by finite quadrupole moment, in 
agreement with mean-field type arguments and the semi-classical theory. In three dimensions, we 
demonstrate that the model in the quadrupolar regime has a phase transition at a finite temperature. 
In contrast to predictions by mean-field theories, the phase transition to the quadrupolar phase turns 
out to be of the second order. We also examine the critical behavior in the two marginal cases with 
the SU(3) symmetry. 

PACS numbers: 75.40.Mg, 75.10.Jm 
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Spin interactions of the order higher than second has 
been discussed for many years iQ p|. There are various 
sources of the high order terms. For example, they may 
arise from the effect of crystalline fields, or the high or- 
der perturbations of electron exchanges. These high or- 
der terms were identified or speculated to be responsi- 
ble for some of the phase transitions observed in various 
magnetic materials ^, |[ ^, |^, ||, ^. In contrast to 
the second order or bilinear interaction models, quantum 
spin models with the high order terms can have a phase 
diagram qualitatively different from their classical coun- 
terparts. In particular, at zero temperature they may 
have non-magnetic ordered phases such as quadrupolar 
phase phase. These non-magnetic phases have been a 
focus of attention in recent years ||^. 

In order for a higher order term to have a non-trivial 
contribution to the Hamiltonian, the spin must be larger 
than or equal to unity. Among the simplest 5 = 1 cases, 
we consider the model with the highest symmetry since 
it probably provides us with a good starting point for 
developing a complete study of wider range of models 
with lower symmetry. In the present article, therefore, we 
discuss 5 = 1 isotropic bilinear-biquadratic Heisenberg 
model: 
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Since the biquadratic term in this Hamiltonian arises 
from the fourth order perturbation of electron exchanges, 
it is usually smaller than the bilinear term that results 
from the second order perturbation. However, it was 
pointed out[^ that the bilinear term can be comparable 
with or smaller than the biquadratic one as a result of 
cancellation of ferromagnetic and antifcrromagnetic con- 
tributions, when we take various hopping terms into ac- 
count. 

For the one-dimensional number of exact solu- 

tions and high-precision numerical calculations have es- 



tablished the character of most of phases and the transi- 
tion points. For the two- or higher-dimensional cases, on 
the other hand, our understanding largely depends upon 
mean- field type approaches or semi-classical theories JTT|. 
A phase transition to non-magnetic ordered phase was 
predicted for a wide range of biquadratic models includ- 
ing the present model. The mean- field approximation [|2| 
was applied to the present model resulting that there is a 
first-order phase transition from the paramagnetic phase 
to the quadrupolar phase (or the spin-nematic phase) 
when the biquadratic interaction is sufficiently large. 

Since the mean-field type approaches are usually ac- 
companied by uncontrollable errors, the confirmation 
through rigorous proof or numerical calculations is neces- 
sary. In the classical model (5 = oo) with small Jl/ Jq, it 
was rigorously proved ||l^ that the quadrupole moment is 
finite in some temperature range above the dipolar transi- 
tion point. In the quantum case of 5 = 1, the quadrupole 
moment was proved [[l^ to be finite at zero-temperature 
in some range of the parameter Jl/Jq in three dimen- 
sions. The range where this rigorous proof applies is 
not the same as, but smaller than the the quadrupolar 
region predicted by the mean-field arguments. In two di- 
mensions, there is no rigorous proof of existence of the 
quadrupolar phase. 

We reported in the previous work[|l5j that the parame- 
ter space of positive Jq is divided by the two SU(3) points 
into three regions; ferromagnetic (— tt < < — 37r/4), 
antifcrromagnetic (— 7r/2 < < 0), and non-magnetic 
(— 37r/4 < 9 < — 7r/2) regions. Here 9 is defined by 



Jl = — J cos 9, Jq 



-Jsm9 



(J>0). 



The nature of the non-magnetic phase was not numeri- 
cally identified in the previous work, although the mean- 
field theory predicted that it is the quadrupolar phase. 
In the present paper, we show for the model with nega- 
tive 9 that (1) the non-magnetic phase is characterized 
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by the finite quadrupole moment in two and three di- 
mensions, (2) a phase transition to quadrupolar phase 
occurs at a finite temperature in three dimensions, and 
(3) the quadrupolar transition is of the second order in 
contrast to the mean-field prediction. We also discuss 
the critical behavior of the three-dimensional system at 
finite temperature. 

In the classical counterpart of the present model, the 
long range order at zero temperature is always dipolar, 
i.e., ferromagnetic or antiferromagnetic, except for the 
special case of 6* = — 7r/2, where the dipolar degrees of 
freedom are non-interacting and disordered. In contrast, 
for the quantum model for S* = 1, it is argued based on a 
mean-field approximation that there is an intermedi- 
ate phase between the anriferromagnetic region and the 
ferromagnetic region, and that this phase is characterized 
by finite quadrupole moment. Because of the limitation 
of the mean-field-type theory, it always predicts a finite 
temperature phase transition to the quadrupolar phase 
regardless of the dimensionality. This is of course wrong 
in one dimension. In two dimensions, too, the existence 
of finite temperature phase transition is very question- 
able because of the Mermin- Wagner theorem. Even at 
zero temperature, the existence of the finite quadrupole 
moment is not totally clear. Mathematically rigorous ar- 
guments, so far, has not established any long-range order 
in the intermediate parameter region. 

In order to answer to the question concerning the exis- 
tence of the quadrupole order, we performed Monte Carlo 
simulation using the loop algorithm proposed in the pre- 
vious letter ijl^. The algorithm removes the ergodicity 
problem and considerably reduces the critical slowing 
down. The energy E, the dipole moment (i.e. magne- 
tization) Mz, the staggered magnetization N^, and the 
quadrupole moment, were measured. We consider only 
the zz component of the quadrupole moment in this ar- 
ticle, which we denote by Qz', 

- E ((SD' - 1) ■ 

The equal-time structure factors and the susceptibilities 
associated to these quantities were also measured. The 
system size ranges from L = 4 up to L = 128 for two- 
dimensional case and up to L = 64 for three-dimensional 
case. For each data point, we typically run the simulation 
for more than 10000 Monte Carlo Steps. 

For each system size in two dimensions, the thermal 
average of the absolute value of the quadrupole moment, 
q = (IQzD/A^ converges to a certain finite value as the 
inverse temperature K = J /k^T increases. Here, the ab- 
solute value \Qz\ is taken in the representation basis in 
which Qz is diagonalized. For any finite system, the con- 
vergence is exponential with some characteristic (imag- 
inary) time scale. Although this characteristic time is 
larger for larger systems, the size dependence is weak. 
Therefore, we can extrapolate the data to the limit of 
K = oo without examining extremely low temperatures. 
After taking the zero-temperature limit numerically, we 
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FIG. 1: The quadrupole moment at zero temperature in two 
dimensions. 

then take the infinite system size limit. The system size 
dependence is algebraic; 

q{L, K — oo) ^ q{L — oo, K — oo) + a/ L. 

This system size dependence is the same as that of the 
staggered magnetization in the antiferromagnetic Heisen- 
berg model in two dimensions. 

Quadrupole moment at zero temperature as a function 
of 9 for various system sizes is plotted in Fig.|l|, together 
with the extrapolation to infinite size. We now see 
that the quadrupole moment is finite in the intermedi- 
ate phase as well as in the dipolar phases. In addition, 
it exhibits discontinuity at the two SU(3) points. Since 
the quadrupole moment is finite whenever the dipole mo- 
ment is finite, it falls down to a finite value, not to zero, 
as we pass the phase boundary from the intermediate re- 
gion to one of the two dipolar regions. Since the dipole 
moment is vanishing in the intermediate phase as we saw 
in the previous paper [Tsl, the quadrupole moment is the 
characterizing order parameter for this phase. 

In order to check the existence or absence of a phase 
transition at a finite temperature, we have examined the 
specific heat. We have observed a broad peak at the 
temperature that roughly corresponds to the saturation 
temperature of the quadrupole moment. The peak height 
and width do not show a significant size dependence, 
indicating that it is not a phase transition but a point 
where a gradual change from the paramagnetic state to 
the quadrupolar state takes place. 

We plot in Fig.^ the size dependence of the quadrupole 
moment as a function of system size, at various temper- 
ature in the case oi 9 — — 0.67r. In this case, the peak in 
the specific heat is located aX K ^ 1.2. In Fig.^j, we see 
that the quadrupole moment shows the asymptotic size 
dependence 

q oc 

down to the temperature K ^ 1.6. For the temperature 
lower than 1.6, the largest system size that we examined 
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FIG. 2: The quadrupole moment as a function of the system 
size ai 6 — ~0.6n in two dimensions. 



is not large enough to see the asymptotic behavior. The 
transition temperature of the KT-type phase transition 
is usuaUy about 10 or 20 percent smaUer than the peak 
temperature of the specific heat. Therefore, if there were 
a Kosterhtz-Thouless type transition, we should be able 
to see a non-trivial algebraic decay for K < 1.6, which 
we do not detect in Fig.^. This indicates that there is no 
phase transition at any finite temperature. 

For the system in three dimensions at zero tempera- 
ture, we again observe three parameter regimes; ferro- 
magnetic, quadrupolar, and antiferromagnetic, with ex- 
actly the same phase boundaries as those in two dimen- 
sions. Namely, the nature of the ground state changes 
at the two SU(3) points, 9 = — 7r/2, —3tt/A. To see this 
in detail, we analyze the order parameters as in the two- 
dimensional case; the extrapolation to zero temperature, 
and then to the infinite system size. The behavior of the 
zero-temperature quadrupole moment as a function of 9 
is similar to the two-dimensional case, but the conver- 
gence to the infinite size limit is faster. The quadrupole 
moment shows discontinuity at the two symmetric points. 
The zero-temperature phase diagram in three dimensions 
turns out to be exactly the same as that in two dimen- 
sions. We speculate that this is true for any dimensions 
except for one dimension. 

Having seen the long range order at zero temperature 
in the intermediate quadrupolar regime, we now ask if 
there is a phase transition at a finite temperature. Even 
in two dimensions, we have seen a broad peak in the spe- 
cific heat and a cross-over behavior from completely dis- 
ordered states to partially ordered states as we decrease 
the temperature. This may be regarded as a precursor 
to the phase transition in higher dimensions. In fact, in 
the specific heat as a function of the temperature, we 
see a much sharper peak in three dimensions than in two 
dimensions. The peak is not only sharp but also shows 
clear size dependence, indicating a phase transition. 

We can clearly see a strong correction to scaling, espe- 
cially in Fig.|[ According to the mean- field theory [|l2j, 
this is a first order phase transition. If this is the case, 
the peak height and width should be proportional to L'^ 
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FIG. 3: Scahng plot of the specific heat at 
three dimensions 
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FIG. 4: Scahng plot of the quadrupole moment at 6 — —O.Gtt 
in three dimensions 



and L"*^, respectively. In other words, a finite size scaling 
plot with exponents d and —d for the vertical scale and 
the horizontal scale, respectively, should work. Other 
quantities should also obey similar scaling forms with 
trivial exponents. We find that this is obviously not the 
case for any quantity. Instead, we assume the following 
finite-size-scaling forms; 

c = L'^y-^c{{K - Kc)Ly), 

q = Li'yq{{K-K,)Ly). 

where c is the specific heat. The best plots are obtained 
with 

y = 1.97 and = 0.6346, 

for the specific heat, and 

y = 2.30, yP = 0.44 and Kc = 0.6345, 

for the quadrupole moment. The scaling plots are 
shown in Fig.|| for the specific heat and in Fig^ for the 
quadrupole moment. The discrepancy among the esti- 
mates of indices may be due to a relatively large contri- 
bution of the non-singular part to the specific heat. We 
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FIG. 5: The phase diagram of the three-dimensional model. 
Curves are mere guides to eyes. 

have estimated the critical temperatures and indices at 
6' = — O.Ttt in a similar fashion, and found that the critical 
indices are close to the corresponding ones for 9 = — O.Gtt 
quoted above. This fact suggests that they belong to 
the same universality class, as expected. Based on these 
results, we conclude 

y = 2.15(20), and y(3 = 0.46(4) 

for -37r/4 < < -7r/2. 

The two SU(3) points, 9 = -7r/2 and -37r/4, are 
of special interest, since the universality class of the 
critical point may be different from the one discussed 
above due to the higher symmetry. For these points of 
higher symmetry, we obtained better scaling plots than 
Fig.^ and Fig.^. The estimated critical temperatures are 
Kc = 0.6389(3) for 9 = -it/2 and = 1.0724(4) for 
9 = — 37r/4. The critical indices for these two cases agree 
with each other, yielding 

2/ = 1.82(5), 2//3 = 0.48(1). 



These result suggest that the critical points of the two 
SU(3) models belong to the same universality class and 
it is distinct from the one for the less symmetric cases 
although the difference in the indices is small. In Fig||, 
we summarize the estimated critical temperatures in the 
form of a 9-T phase diagram. 

To summarize, we have studied the isotropic bi- 
quadratic Heisenberg model in two and three dimensions 
for negative 9. In two dimensions, we have identified 
the intermediate phase as the quadrupolar phase. The 
phase transition at a finite temperature has been ex- 
cluded. In three dimensions, we have studied finite tem- 
perature properties as well as zero temperature ones. At 
zero temperature, the phase diagram is exactly the same 
as that of the two-dimensional case. We have found that 
there is a finite temperature phase transition not only in 
the ferro- and antiferro- magnetic regimes but also in the 
quadrupolar regime. In contrast to the mean-field predic- 
tion, the transition to the quadrupolar phase has turned 
out to be of the second order. The critical indices are 
also estimated. While the two SU(3) symmetric points 
belong to the same universality class, it is suggested to be 
distinct from the one for the less symmetric (i.e., SU(2)) 
models. Studies on the properties of low- lying excitations 
are still in progress and will be reported elsewhere. Less 
symmetric models with higher order interactions may be 
more important than the present model from the prac- 
tical point of view, since the higher order interactions 
in real magnets often arise from the crystalline effects, 
which have lower symmetry. Studies on some of these 
models are also in progress. 
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